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We have investigated the effects of inplane and interplane nearest neighbour overlap integrals (so 
and si) and site energy difference between atoms in two different sublattices in the same graphene 
layer (A) on the electronic dispersion of bilayer graphene within tight binding model. Also, mod- 
ifications in bilayer graphene bands due to inplane next nearest neighbour interactions (-yu, Su) 
and next to next nearest neighbour interactions (721, S2i) have been studied. It is observed that si 
introduces further asymmetry in energy positions of top conduction band and bottom valence band 
at the K point on top of the asymmetry due to A. Moreover, so, Si as well as the other inplane 
coupling parameters induce noticable electron-hole asymmetry in the slope of the bands and changes 
in band widths. The density of states of bilayer graphene has also been calculated within the same 
model. 

I. INTRODUCTION 

Bilayer graphene is a system of two stacked hexagonal graphene sheets. Among all the carbon based materials 
of recent interest, bilayer graphene is of very much importance because this is the only two dimensional material 
in which the band gap between valence and conduction bands can be controlled by applying an external electric 
field perpendicular to the laycrs^^^ or by chemical doping of one of the layers^. This makes it a potential 
candidate for future application in nanoclcctronics. There have been an intensive research work going on over the 
last several years to understand the intresting behaviours it manifests. Since electronic band structure is one of the 
key ingredients to shed light on the understanding of the properties of a material, in particular, as far as electronic 
properties of the material are concerned, a significant amount of effort has been put to study the band structure of 
bilayer graphene&i^^. Regarding this the major concern was to see how the linear dispersion near the K ponit in 
single layer graphene is modified in presence of interlayer coupling in bilayer graphene. It is now well established that 
even the slightest presence of the interlayer coupling kills the linearity of the dispersion of monolayer graphene and 
converts it to a parabolic one for the bilayer graphene. Moreover, most of the studies have talked about the symmetric 
nature of the valence and conduction bands of bilayer graphene. Only some very recent work o 10 ' 11 ! 12 have reported 
about asymmetry between the valence and conduction bands in this system. In this work we have investigated for 
some other parameters which also contribute appreciably towards asymmetry in bands of bilayer graphene. While the 
existing literature mostly discussed about interlayer coupling energies in bands of bilayer graphene and in graphite^, 
we here report about asymmetry in band structure of bilayer graphene due to the inclusion of inplane and interplane 
nearest neighbour overlap integrals (sq and s\) in the tight binding band structure calculation because presence of 
nearest neighbour overlap integral (so) gives quite a bit asymmetry in the band structure of monolayer graphen o ' 15 . 
We have also studied the effect of inplane next nearest neighbour (nnn) and next to next nearest neighbour (nnnn) 
hopping energies (7^ and 72^) along with the corresponding overlap integrals (si, and S2i) on the band structure of 
this system in line with a similar work on single layer graphene^. To calculate the bands we have used A = 18 meV 
and nearest neighbour interplane coupling energy 71 = 400 mcV— . For all the inplane parameters, we have used a 
set of parameters which we have determined by comparing the tight binding bands of single layer graphene with first 
principle results by considering the fact that magnitudes of the parameters should decrease with increasing distance 
in our earlier workii. The article has been arranged in the following way: in section I we talk about the structure 
of bilayer graphene; section II focusscs entirely on electronic structure of bilayer graphene with nearest neighbour 
inplane and interplane couplings, inplane next nearest and next to next nearest neighbour couplings; section III shows 
the density of states of bilayer graphene for the above mentioned bands and finally, summary and conclusions over 
the previous sections are given in section IV. 



A. Structure 



Bilayer graphene is a coupled system of two monolayers of graphene with usually an AB-stacking fashion of the 
layers, i. e. the atoms in one layer arc not just on top of the corresponding atoms in the other layer, rather the 
arrangement is such that if one layer is projected on the other, the A-type (say) atoms coincide with the A-type 
atoms but the B-type atoms come at the center of the hexagons of the other layer—. So A and B type atoms should 
not be treated as chemically equivalent carbon atoms because they belong to two different chemical environments. 
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But since this site energy difference is very small (~ 18 meV), they are frequently treated as chemically equivalent 
carbon atoms. The structure of AB stacked bilayer graphene and its brillouin zone are shown in figure (1). The 
distance between two atoms in a layer is 1.42A and the interlayer separation is 3.35A. The brillouin zone of bilayer 
graphene is same as that of single layer graphene because it does not have periodicity along z-dircction. Bilayer 
graphene unit cell contains four atoms. The relevant hopping parameters, i. e. nearest neighbour inplane coupling 
energy (70), nearest neighbour interplane coupling energy (71), next nearest neighbour inplane coupling energy (711) 
and next to next nearest neighbour inplane coupling energy (721) have also been shown in the figure. 
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FIG. 1: (a) Structure of bilayer graphene with unit vectors Si, S2; the intralayer nn coupling energy (70), interlayer nn coupling 
energy (71), intralayer nnn coupling energy (71;) and nnnn coupling energy (721)- 00 Brillouin zone of bilayer graphene with 
unit vectors 61, 62 and the high symmetry directions. 



II. ELECTRONIC DISPERSION OF BILAYER GRAPHENE 



In the following first we will be talking about some general aspects of the formalism of tight binding band structure 
of bilayer graphene, then the details of tight binding band dispersion under different situations like: (i) the simplest 
case with nearest neighbour inplane and interplane hopping of ir electrons; (ii) modifications in the bands due to 
the inclusion of overlap integrals coming from the same neighbours and due to site energy difference between A and 
B atoms in two different sublattices in the same layer; (ii) effect of inplane second nearest neighbour hopping and 
corresponding overlap integral; (iv) effect of inplane third nearest neighbour hopping and overlap integral on the 
band near the brillouin zone corner (K) point and over the whole brillouin zone, and the density of states of bilayer 
graphene corresponding to the above mentioned dispersions. Since bilayer graphene unit cell contains four atoms 
coming from two different sublattices of two layers, the Bloch wave functions for A and B type atoms for each layer- 
are 



tf£*(r) = 1/VN ®A(r - r Ai )e ik ' rAi 



Ai=i- 



ttf (r) = 1/VN J2 <$>B(r-r Bi )e lk - r ^, 



where i refers to layer index, $ s are p z atomic orbitals, N is the number of unit cells in the crystal, Ca { and Cb { 
are contributions to the total wavefunction coming from i th atom of A and B sublattices respectively. The total wave 
function is 



If H is the Hamiltonian, the general secular equation including overlap integrals becomes 



Ha 1 A 1 ~ E(k)SA 1 A 1 H Al Bx - E(k)SAxBx H Al A 2 - E(k)SAxA 2 

HbxAx - E(k)S Bl Ax H Bl Bx - E(k)S Bl Bx H Bi a 2 - E(k)S Bl A 2 

H-A 2 Ax — E(k)SA 2 A x El A2 Bx — E(k)SA 2 Bx El A2 A 2 — E(k)SA 2 A 2 

Hb 2 Ax — E(k)SB 2 A 1 Hb 2 Bx — E(k)SB 2 B! Hb 2 a 2 — E(k)SB 2 A 2 



H Ai b 2 - E(k)S Al B 2 

HbxB 2 — E(k)SBxB 2 

Ha 2 b 2 — E(k)SA 2 B 2 
Hb 2 b 2 - E(k)SB 2 B 2 



0, 



(1) 
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where i?„ s „. and S r 



are defined as 



S x , 



tij 



(*fc*l*fc*) respectively. 



Here x, y stand for both A and B; i, j can take either of the values 1 and 2, and S s are the overlap integrals. 



A. Modifications in the bands due to the inclusion of overlap integrals coming from inplane and interplane 
nearest neighbours and due to breaking of sublattice symmetry 

In this part, the band structure of bilayer graphene has been calculated in presence of nearest neighbour inplane 
and interplane coupling energies (70, 71) along with the corrections coming from the corresponding overlap integrals 
(s , si). The sublattice energy difference (A) has also been included. When solved, the eigenvalue equation readily 
gives the following solutions in case of the simplest possible situation with A, sq and s± to be zero. 



E 1 (k) = -1/2 



7i - v // 7i 2 +47off(fc) 



E 2 (k) = 1/2 



E 3 (k) 



-1/2 



7i + \/7i 2 + 4 7o5(fc) 

7i + \fil+Hg(k) 



and 



Ei{k) = 1/2 



7i - sfli +47o5(fc) 



where g (k) = ^1 + 4 cos 2 (k y a/2) + 4 cos {\/?>k x a/2) cos (k y a/2)j . 
Expansion of the dispersions around the K (27r/v / 3a, 27r/3a) point give 



Ei(k) = -7i-(3/4)(7 2 / 7l )a 2 (Afc) 2 , 

E 2 (k) = -(3/4)( 7 2 /7i)a 2 (Afc) 2 , 

E 3 (k) = +(3/4)( 7 2 / 7l )a 2 (Afc) 2 and 

E 4 (k) = 7 1 + (3/4)(7 2 / 7l ) a 2 (Afc) 2 , 

where Afc is a small change in k around the K point. From the above results we see that there are four it bands in 
bilayer graphene. All the bands are parabolic near the brillouin zone corner. Two of the bands are degenerate with 
zero energy at the K point. -71 and 71 are the energies of the other two bands at that point, i. e. ±71 is the energy 
separation between the highest conduction band and the lowest valence band at that point. Also the valence and 
conduction bands are symmetric over the entire brillouin zone. Since each atomic site has one p z electron, the fcrmi 
energy (Ep) for undoped bilayer graphene is at zero energy. The bands are plotted in figure (2) along with the bands 
with nearest neighbour inplane overlap integral (so), nearest neighbour interplane overlap integral (si) and sublattice 
asymmetric energy (A). Including finite values of A, sq and s± the secular equation becomes 



E a + A-E(k) a-bE{k) 71 - Si-E(fc) 

c-dE(k) E -E(k) 

7l - Sl E(k) E + A - E{k) c - dE{k) 

a- bE(k) E a - E(k) 



0, 



(2) 



where a = 7o/(fc), b = sof{k), c = 70/* (k) and d = sof*(k). Following are the eigen solutions to the above equation. 



-B 1 ±,/S 2 -4A 1 C 1 



where A 1 = 1-bd+sx, B 1 = ad + be - 71 - s 1 E Q - 2E Q - A, C\ = jtE - ac+ E (E + A) 
and£ 3 ,4( fc ) = -B 2 ± J B% - 4A 2 C 2 /2A 2 , 



(3) 
(4) 



where A 2 



= l-bd-sx, B 2 = ad + be + 71 + s 1 E Q - 2E Q - A, C 2 = 



-71 £0 -ac + E (E + A). 
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FIG. 2: Electronic spectra of bilayer graphene. For all three cases, in the upper panel we have plotted the bands along all 
the high symmetry directions of the brillouin zone and the lower panels show zoomed versions of the corresponding upper 
panels near the K point. In all the panels the black curves represent the symmetric bands, i. e. bands in presence of nearest 
neighbour inplane (70) and interplane (71) hopping energies. In (a) the green curves contain informations due to the presence 
of sublattice asymmetry (A), clearly visible from the lower panel; in (b) the blue curves have the effect of nearest neighbour 
inplane (so) and interplane (si) overlap integrals but zero A, whereas the red curves are plotted with non-zero values of so, 
si and A. Effects of these parameters at the K point are visible from the lower panel. In (c) the bands in presence of finite 
values of so, si and A (red curve) are compared with the symmetric bands (black curve) over the whole brillouin zone (upper 
panel) and nearer to the K point within optical energy range (lower panel). 

We first discuss the effect of the asymmetry term (A) on the spectra at K point. The two bands which are having 
degeneracy at K point are not affected by A at that point but the other two bands appear at energies A — 71 and 
A + 71, i. e. though the separation is still ±71, they become asymmetric regarding energy positions at that point. This 
term has negligible effect on the slope of the curves. The lower panel of (a) in fig. (2) shows that the lower valence 
band comes closer to Ep whereas the upper conduction band moves away from Ep when the sublattice asymmetry is 
taken care of. Next we illustrate the effect of overlap integrals sq and si at K point. It is observed that the inplane 
nearest neighbour overlap integral so has no effect at K point because it always appear in product with g(k) which 
vanishes at K. But the interlayer overlap integral si, though not altering the degenerate bands, affects those bands 
lying further from Ep quite significantly. The upper lying band shifts further to —71/(1 — s±) and the lower lying 
band comes closer to 7i/(l + si) (blue curve in lower panel of (b) in fig. (2)). The energy separation between them 
becomes 271/(1 — sf). When A and overlap integrals are considered together, the upper lying band shifts even further 
to (A — 7i)/(l — Si) and the lower lying band comes more close to (A + 7i)/(l + si) (red curve in lower panel of (b) in 
fig. (2)). Under this condition the above energy separation becomes 2(71 — siA)/(l — sf). Moreover, sq and Si play 
important role in changing the slope of the curves. They push off the conduction bands and pull in the valence bands 
almost over the entire energy range, particularly near the T point the effect is most prominent (red curves in upper 
panel of (c) in fig. (2)). With a close look over a smaller energy range (lower panel of (c)) we see that the conduction 
bands are repelling each other towards higher energy side whereas the valence bands have the tendency to merge. 

B. Effect of inplane second nearest neighbour hopping and corresponding overlap integral on the spectrum 

of bilayer graphene 

In this part we have studied the changes in the above dispersions by taking care of inplane second nearest neighbour 
transfer integral (71^) and the corresponding overlap integral (su). Since the inplane next nearest neighbour atoms 
belong to same sublattice in a plane, only the entities Ha 1 a 1 , Sa 1 a 1 ] EIa 2 a 2 , Sa 2 a 2 ] Hb 1 b 1 , Sb 1 b 1 and Hb 2 b 2 , Sb 2 b 2 
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get modified and changes occur in A\, Bi, C\ and A2, B2, C2 as follows 

At = (1 + s u u{k)) 2 - bd + si(l + s u u{k)), 

B x = ad + bc- 71(1 + s u u(k)) - s^Eq + 7 H u(fc)) - (1 + s u u{k))(2E + 2>y u u{k) + A), 

Ci = i 1 (E +i u u{k))-ac+{E +i li u(k))(E +'y li u(k) + A), and 

A 2 = (l + s u u(k)) 2 -bd- Sl (l + s liU (k)), 

B-2 = ad + bc + 7l (l + s u u(k)) + s^Eq + lu u{k)) - (1 + s u u(k)){2E + 2 lu u{k) + A), 

C 2 = -ji(E +juu(k))-ac+(E + 'Yuu(k)){E +juu(k) + A), 

where u(k) = 2 cos (k y a) + 4 cos (k x aV3) cos (k y a/2). 

C. Modification due to inplane third nearest neighbour hopping energy and overlap integral on the bands 

of bilayer graphene 

Inplane third nearest neighbour atoms belonging to the other sublattice incorporate changes to the quantities 
Ha 1 b 1 i Sa 1 b 1 ', Ha 2 b 2 i Sa 2 b 2 an d their complex conjugates which inturn modify A±, B±, C\ and A2, B2, C2 through 
changes in a, b, c and d. Third nearest neighbour inplane transfer energy and overlap integrals are denoted as j 2 i and 
S2i- Modified a, b, c and d appear as 

a = 7o/(&) +l2ifi(k), b = s f(k) + s 2l fi(k), c = 70/* (k) + l2ifZ (k) and d = s Q f*(k) + s 2i fi(k), 
where = e ifc Wv^ 2 cos k y a + e - 2lk * a ^. 

So the eigen solutions in this case are same as equation (3) and (4) with a different set of A\, B\, C\ and A2, B2, 
C2 due to differences appearing in a, b, c and d. To illustrate the effects of inplane next nearest neighbours and next 
to next nearest neighbours on the spectra of bilayer graphene, all the bands are plotted together with the symmetric 
spectra (black curves) and the spectra with corrections due to nearest neighbour inplane and interplane overlaps (red 
curve) in figure (3). From the plots shown in fig. (3) and from plots not shown here it is very clear that presence 
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FIG. 3: Electronic dispersion bilayer graphene in presence of nearest neighbour inplane and interplane transfer integrals (black 
curve), nearest neighbour inplane and interplane transfer integrals, overlap integrals and sublattice asymmetric energy (red 
curve), inplane next nearest neighbour interactions (green curve) and inplane next to next nearest neighbour interactions (blue 
curve). The parameters used for these bands are given in Table I. 

of any of the parameters (so, si, 7h, Su, 724, S2i) or all of them together introduces asymmetry of different amounts 
on the symmetric bands due to 70 and 71. While A makes the upper conduction and lower valence band positions 
asymmetric at K point, others give asymmetry to the bands almost over the whole brillouin zone. We have already 
discussed the effects of so, s± on the spectra with 70 and 71 in the previous section. Now, we describe the effects due 
to inplane next nearest neighbours (71;, Sh) on top of sq, si. There is no significant change in energy positions of the 
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top conduction band and the bottom valence band at the K ponit in presence of inplane nnn and nnnn interactions. 
For the valence bands the nnn bands are nearer to the symmetric bands within ~ 2 eV but the nnnn bands are closer 
to the modified bands due to Sq and s\. For the conduction bands the scenario is opposite within the same energy 
range. Beyond that the nnn bands move away from the symmetric bands, though in opposite directions for valence 
and conduction bands with respect to fermi energy, nnnn bands start moving towards the symmetric bands such that 
the nnn and nnnn valence bands almost coincide and the co duction bands keep a small separation among themselves 
around T point . The values of the parameters used for the plotted bands are given in Table I. In presence of inplane 



TABLE I: Tight binding parameters 



Curves 


£o(eV) 


A(eV) 


70 (eV) 


7i (eV) 


7ii(eV) 


72i(eV) 


so 


si 


Sli 


82i 


Black 






-2.78 


-0.4 














Red 




0.018 


-2.78 


-0.4 






0.117 


0.04 






Green 


-0.45 


0.018 


-2.78 


-0.4 


-0.15 




0.117 


0.04 


0.004 




Blue 


-0.45 


0.018 


-2.78 


-0.4 


-0.15 


-0.095 


0.117 


0.04 


0.004 


0.002 



next nearest neighbour interactions there arc three choices of Eq (= i^u 1 i^u — A ± 71) for which the fermi energy 
comes at zero. For other choices of Eq the the fermi level will shift from zero. In presence of inplane next to next 
nearest neighbours the choices of E for which fermi energy will be at zero are same as in case of nnn interactions. 

D. Density of states of bilayer graphene 

The density of states (DOS) of bilayer graphene corresponding to the bands in fig. (3) with the parameters in Table 
I are shown in figure (4). From the curves it is evident that the finite, though small, density of states (DOS) at 




FIG. 4: Density of states of bilayer graphene in presence of nearest neighbour inplane and interplane transfer integrals (black 
curve), nearest neighbour inplane and interplane transfer integrals, overlap integrals and sublattice asymmetric energy (red 
curve), inplane next nearest neighbour interactions (green curve) and inplane next to next nearest neighbours (blue curve). 
The parameters used in these curves are given in Table I. 

the fermi energy of bilayer graphene remain unaffected due the presence of factors leading to asymmetry in bands. 
Further, the slope of the density of states curves within < 0.5 eV does not change due to the above mentioned factors. 
Rather, they have prominent effect in bringing asymmetry in band widths. To illustrate, we sec that when only 
nearest neighbour inplane and interplane transfer energies (70 and 71) are there, valence band DOS and conduction 
band DOS (black curve) are exactly symmetric with respect to the fermi energy at zero, i.e. band widths arc same. 
Also the van-Hov singularities are at symmetric positions, the scenario being consistent with the corresponding band 
structure. As soon as the nearest neighbour inplane and interplane overlap integrals (so and s\) are taken into account 
the valence band (VB) DOS and conduction band (CB) DOS start becoming asymmetric, i.e. the band widths become 
different: valence band becomes narrow and conduction band widens. Moreover, the van-Hov singularities in VB come 
closer to Ep whereas those in CB move away from Ep. In presence of inplane second nearest neighbours CB becomes 
slightly narrow and VB slightly wide compared to the previous case, though VB is still narrower compared to CB. 
There is no significant change in the positions of the van-Hov singularities in this case when compared with the 
previous one. Inplane third nearest neighbours do not have much effect on the widths of the bands on top of that of 
the inplane second nearest neighbours but bring the van-Hov singularities slightly nearer to the fermi energy. It is 
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also observed that the sublattice asymmetric energy (A) does not have any significant effect on the density of states 
of bilayer graphene. 

III. SUMMARY AND CONCLUSIONS 

To summarize, we have illustrated the effects of various parameters governing electron-hole asymmetry in the band 
structure of bilayer graphene within tight binding model. When compared the role of site energy difference (A) 
between A and B sublatticcs in the same graphene layer on the electronic spectra of single layer graphene and bilayer 
graphene, a distinct difference is observed between the two systems. Sublattice asymmetry in monolayer graphene 
introduces a gap in the spectra at the K point whereas in bilayer it does not induce gap in the spectra, rather it gives 
an asymmetry in the energy positions of the top valence and bottom conduction bands with respect to the energy 
at which the other two bands are degenerate. Moreover, in presence of A the gap between top valence and bottom 
conduction bands at the K point remain intact to ±71 which is the separation even without A. Apart from A the 
other important factor which contributes significantly to the positions of the top valence and bottom conduction bands 
at K point is the interlayer nearest neighbour overlap integral (si). Hence, we find that E c (top) and ^(bottom) are 
functions of A and si both at K point. Regarding this ref. (11) discusses only about the dependence on A. Our 
study suggests for the consideration of a model containing both A and si for a more accurate determination of A 
from experimental results. Further, we observe a considerable change in the slope of the bands in presence of nearest 
neighbour inplane and intcrplanc overlap integrals (sq and s\) compared to those with nearest neighbour inplane and 
interlayer coupling energies (70 and 71) only. Ref. (11) has discussed about the induction of electron-hole asymmetry 
in the slope of valence and conduction bands due to next nearest neighbour interlayer coupling energy (74) but without 
overlap integrals. It is noted there that with finite 74 two conduction bands are closer and the valence bands are 
further apart at k values away from the K point compared to the bands without 74. We observe that finite values of 
so and si also give electron- hole asymmetry in the slope of valence and conduction bands but with an opposite trend 
compared to the bands with 74. With finite values of sq and Si the conduction bands move futher from each other 
and the valence bands come closer for momentum values away from the K point value. It could be concluded from 
this comparative study that the combined effect of 74, sq and si could be somewhat balancing between the above 
two cases, may not be a complete balance to get back the symmetric spectra with only 70 and 71 but the degree 
of electron- hole asymmetry in slope of the bands will get modified. A model including 74, Sq and si, though very 
complicated to handle, may lead to more accurate determination of the important parameter like 71 when comparing 
experimental results having asymmetry in electron and hole sides. Now comes the inplane next nearest neighbour 
interactions (71.; and su). With and su the trend of electron-hole asymmetry in slope of bands remain similar as 
that with so and si but the asymmetry in valence and conduction band widths near T point is reduced compared to 
that with zero values of ju and Su- Moreover, in this case if the site energy term Eq is not chosen properly there 
is a shift in fermi energy^. Inplane third nearest neighbour interactions (721 and S2i) do not affect much on top of 
inplane second nearest neighbours except very little change in band withs at T point and slight modifications in the 
slope of the bands. Hence, as far as electron-hole asymmetry in slope of valence and conduction bands are concerned, 
a Hamiltonian including so, si and inplane second nearest neighbour interactions (ju and Su) could be sufficient to 
interpret experimental results (e. g. cyclotron resonance data) with asymmetry in electron and hole sides. Also, the 
inplane third nearest neighbour interactions could be more useful in determining all the above mentioned parameters 
by fitting bilayer graphene bands with first principle results or with angle resolved photoemission(ARPES) data over 
the whole brillouin zone. 
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